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ABSTRACT 

Inclusion of the inductive electric field, Eind due to the temporally changing B, in magnetic explo- 
sions is discussed, with emphasis on solar flares. Several roles played by Eind are identified: on a global 
scale, Eind produces the EMF that drives the explosion; the associated Eind x B drift is identified with 
the inflow of magnetic field lines into a reconnection region; the polarization current, associated with 
dEind/dt, implies a J x B force that accelerates this inflow; and the component of Eind parallel to 
B accelerates the energetic electrons that cause hard X-ray emission and type III radio bursts. Some 
simple models that describe these effects are presented. A resolution of the long-standing "number 
problem" in solar flares is suggested. 

Subject headings: Sun: flares; Magnetic fields; Magnetic reconnection; Acceleration of particles 



1. INTRODUCTION 

As discussed in an accompanying paper, referred to as paper I, magnetic explosions are intrinsically time-dependent 
phenomena, and must involve an inductive electric field, Eind, and a displacement current Jdisp = ^od^inddt. These 
fields are not taken into account explicitly in conventional models for solar flares and other magnetic explosions; the 
physics associated with them is either treated in terms of quasi-electrostatic proxies or ignored. The following six roles 
for Eind and Jdisp need to be taken into account in models for magnetic explosions. 

1. On a global scale, the electromotive force (EMF), $ = —d'^B/dt, that drives the explosion-related current is 
determined by the rate of change of the magnetic flux, 5* b , enclosed by the circuit around which this current 
flows. 

2. The electric drift, Uind = Eind x B/|Bp, associated with the inductive field transports magnetic energy into the 
energy-release (or reconnection) region: no mechanical driver is required. 

3. The collective response of the plasma particles to the displacement current results in the polarization current, 
Jpoi = {c^ /v\) Jdisp ±, where va is the Alfven speed. 

4. The J X B force due to the polarization current accelerates Uind, providing the link that allows the changing 
magnetic field to drive the magnetic energy release. 

5. The Alfvenic energy flux that transports energy from the energy-release region to the acceleration region is set 
up by the polarization current changing the current profile in the linking flux tube, and driving a vortex motion 
of the plasma within the flux tube. 

6. Acceleration by E\\ relies on Eind and Jdisp (jSong fc Lvsakl[2006[ ) . 

The first of these roles is discussed in paper 1, and the other five are discussed in this paper. 

A l ong-standing "number problem" concerns the acceleration of the electrons that produce hard X-rays (jMacKinnon fc: Brownl 
Il989| ): the number of accelerated electrons inferred from X-ray data exceeds the number available by a large factor, 
M = 10^ assumed here. A way of resolving the number problem, by including a return current in the model for 
Alfvenic energy transport, is suggested in SJS] 

The argument for the inductive field being the flare driver is presented in 5j2j The role of the polarization current in 
driving this inflow is discussed in fJHl The role of time-dependent fields in Alfvenic energy transport to the acceleration 
region is discussed in 21 Acceleration of electrons by , including the number problem, is discussed in ^ The 
results are discussed and summarized in SjG] 

2. INDUCTIVE FIELD AS A FLARE DRIVER 

In this section it is assumed that the parallel component of the inductive field is screened by charges. The perpen- 
dicular component is interpreted in terms of the E x B drift velocity, Uind- It is argued that Uind can be identified 
as the inflow velocity, caused by the time-changing magnetic field, into the reconnection region. This contrasts with 
a conventional steady-state model is which the inflow is imposed as a boundary condition, and is attributed to a 
mechanical driver. 
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Fig. 1. — A cartoon of the Sweet-Parker reconnection model: the magnetic field (solid) and flow (dashed) lines are confined to the x-z 
plane; dissipation is confined to a current sheet, shown in gray, about the magnetic neutral plane, x = 0. 

2.1. Inductive drift 

The assumption that the parallel component of the inductive electric field is screened implies the presence of a charge 
density, p. This charge density is such that it produces a parallel electric field that is equal and opposite to the parallel 
component, Ei„(j||, of the unscreened Ei„ci- The remaining (perpendicular) component of the inductively-induced 
electric field has a nonzero divergence 



curlEi„d± = 

The inductive drift is the electric drift due to Eind ± '■ 

Eind_L = — Uind X B 

Combining ([!]) and © gives 



divEi„d_L = 



Uind 
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= curl(uind X B), 
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which is the equation for magnetic field lines moving with the velocity Uind- The screening charge density is 

p = eo [uind • curl B - B • curl Uind] ■ (4) 



2.2. Sweet-Parker-type model 

Consider any steady-state reconnection model, such as the Sweet-Parker-type model illustrated in Fig. [1] Such 
a model involves an inflow of plasma and magnetic field into a current sheet, with this inflow being specified as a 
boundary condition. Physically, the infiow must be driven in order for fast reconnection to proceed, and in a st eady- 
state model this requi res a mechanical driver acting across field lines, usually assumed to be a pressure gradient (jSatol 
119781: lHaerendeliri994D . However, the pressure gradient is much smaller that the available Maxwell stre ss, due to the 
J X B force , and the latter is a more plausible driver, as proposed in Uchida's quadrupolar model (jUchidal [19961 : 
iHirose et al.ll200"ll ) for an erupting filament. Here it is suggested that the J x B force acts as a driver under more 
general conditions, provided that the (artificial) steady-state assumption is abandoned. 

A simple model that illustrates the role of Ujnd plays in reconnection involves a magnetic field, that varies globally 
as a function of time, f{t). Although such an assumption is also artificial, it complements a steady-state model. The 
function g{t) = df(t)/dt can be regarded as simulating the rate of reduction in stored magnetic energy beyond the 
boundary of the steady-state region; the associated Uind becomes the inflow velocity. 

In the simplest example of a planar geometry, let the magnetic fleld along the z axis, and oppositely directed above 
and below the current sheet at x = 0. Assuming that the x- and ^-dependences are independent, one may write 
Bz{x,t) = f{t)b{x), with b{—x) = —b{x). Assuming that the system is in uniform in y, one has 

Eindy{x,t) = -g{t) J dx'b{x'), f ^ 

The inductive drift, 

Uindx{x,t) = \ [ dx'b{x'), (6) 

f{t)b{x) Jo 

is along the x-axis, with Uindi:(~^c, t) — —Uindx{x,t). For a decreasing field, g(t) < 0, the flow is directed towards the 
current sheet at x = 0. A generalization to allow the magnetic field lines to be curved in the x-z plane is needed to 
allow outflow along the z axis. Assuming Bx{x,z,t) = f(t)bx{x,z), Bz(x,z,t) = f{t)bz{x, z), the inductive electric 
field is again along the y axis, 

Eindyix,z,t) = -g(t) 



J 

Jo 



dx'bz{x',z)— / dx'bx{x,z') 



(7) 
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and Uind has a component along the z-axis, such that Ujnd • B = 0. In both these simple models there is no parallel 
component of the inductive electric field, so that screening of the parallel component of the inductive field is not 
relevant, and one has p = 0. 

2.3. Self-similar model 

A generalization of the foregoing discussion is possible in the special case where the change in the magnetic field is 
self-similar. Specifically, a self-similar model corresponds to B(x, i) — /(i)b(x). The inductive electric field can then 
be written Eind(x.t) = g(i)e(x). The inductive drift is 

uU^,t)^— i^^^^^i^ . (8) 

The assumed screening of the parallel component of Ejnd implies e(x) • b(x) = and requires p(x, = g{t)r{x.) ^ 0. 
Introducing the vector potential, A(x, t) = g{t)a{x), and the scalar potential, 0(x, t) — g{t)p{x) in the Coulomb gauge, 
diva(x) = 0, one requires that p(x) depend on the distance, s, along field lines, with d^p^x) / ds^ = r(x)/eo- One then 
has 

e(x) = — a(x), b(x) — curla(x). (9) 
The parallel-screening condition becomes 

a(x) • curla(x) 0. (10) 

The flow ([5]) is along the direction of 

a(x) X curla(x) = grad (i|a(x)p) - [a(x) ■grad]a(x). (11) 

The condition (jlOp restricts the form of a magnetic field that can evolve in a self-similar manner. Seeking a solution 
that is separable in Cartesian coordinates, one finds that (ITUl) is satisfied by 

b(x) = {X[Y' - Z'],Y[Z' - X% Z[X' - r']), (12) 

where X — X{x), Y = Y{y), Z — Z{z) are arbitrary functions, and X'.Y'^Z' are their derivatives. The model 
corresponding to ([7]) is reproduced by y = 0, and the model corresponding to ([5]) is reproduced by y = 0, Z = 1. 

3. ROLE OF THE POLARIZATION CURRENT 
In this section, the role of the polarization current associated with a changing inductive electric field is discussed. 

3.1. Polarization current 

The presence of a time-dependent electric field in a pla sma causes a ch arged particle, with charge q and mass m, to 
drift across the magnetic field. This polarization drift is (jNorthrod 1 1 9631 ) 

''P"'-^^- ^^^^ 

In the following, dE^/dt is identified as &Eii^d±ldt, whose curl is determined by d^'B/dt^. The drift (|T3)) depends 
on the sign of the charge, and summing over all charges it implies a current that is proportional to the displacement 
current. This current is referred to as the polarization current: 

Jpoi^^eo^. (14) 

This polarization current may be interpreted as the response of the plasma to the displacement current, and an 
alternative derivation of (jl4p follows by considering the low-frequency response of a plasma. The susceptibility tensor 
for a magnetized plasma at very low frequency is approximately diagonal, with components c^/w^, c^/w^, — Wp/w^, 
where ujp is the plasma frequency. This tensor relates the polarization, P/eqi to the electric field, with the current 
density related to the polarization by J = dP /dt. Differentiating the perpendicular components with respect to time, 
assuming for simplicity that v\ is a constant, and integrating the parallel component (with —uj^ — >■ /dt^) implies 
that the response can be rewritten in the form 

aE^ ajii „ , , 

£o^, -^^eQUpE\\. (15) 



V 



This provides an alternative derivation of Jpoi = J±. This alternative formalism allows one to include dissipation, 
through the conductivity tensor, with perpendicular component and parallel component , leading to additional 
terms, (J_lE_l and a\\ dE\\/dt, on the right hand sides of the two equations (fT5|) . respectively. 
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On differentiating the expression ([2]) for Uind with respect to time, neglecting the time-dependence of B, the result 
becomes 

= Jpoi X B, (16) 

where 77 is the mass density. It follows that the polarization current provides the force that accelerates Ui„d- 

These results lead to the following interpretation of the energy inflow into the energy-release region. Magnetic 
energy is stored in a global current system, and some trigger causes B to start decreasing in time. The associated 
Jpoi accelerates an inductive flow Uind that transports energy over large distances to the (ill-defined) boundary of the 
energy-release region. Within the energy-release region the energy fiow may be described by a steady-state model, 
with further acceleration by Jpoi being unimportant. The rate of decrease in the globally stored energy is balanced by 
the outflow of energy from the energy-release region. 

The conceptual change from a conventional model is that no pressure gradient or other mechanical force is need to 
drive the inflow. As the energy release begins, Jpoi becomes nonzero, and the inflow is driven by the Jpoi x B force. 

3.2. Energy-release region 

The conversion of magnetic energy into the energy in nonthermal particles in a magnetic explosion is assumed to 
occur in two stages: conversion of magnetic energy into an Alfvenic flux in an energy-release region, and conversion of 
the Alfvenic flux into energetic particles in a remotely located acceleration region. The energy conversion in the first 
of these region is discussed here. 

Conversion of magnetic energy into another form requires some dissip ation, but the dissipati on need not be energet- 
ically significant. Reconnection is effective only around a magnetic null (jPriest fc Forbesll200ClD . and the Sweet-Parker 
model illustrated in Fig. [T] shows the inflow and outflow in such a region. To be effective in a solar flare, reconnection 
must be fast (compared with the Sweet-Parker rate), and the geometry is similar for fast reconnection models. Numeri- 
cal simulations show t hat fast (Petsche k) reconnection can o ccur only if the resistivity is localized (Biskamo & Schwar3 
[19971 lUzde nskv"2003') and anomalous. (jKHEmSIMl EOO^, in which case the reconnection region is on a micro scale 
(pape r 1). A statistically la rge number of such reconnection sites is needed, and this requires multiple local magnetic 
nulls (iLazarian fc Vishniad lT999) . Although dissipation is an essential ingredient in reconnection, the energy dissipated 
in converting the magnetic energy into kinetic energy or an Alfvenic flux is arbitrarily small for arbitrarily small regions 
of anomalous resistivity. It is assumed here that magnetic energy conversion occurs without significant dissipation in 
the energy-release region. 

The energy flow into the energy-release region is inductively driven, as discussed above. Within the energy-release 
region a steady-state model relates the energy outflow to the inflow. The outflow can be in three forms: kinetic energy, 
magnetic energy transported due to fleld lines frozen into the flow, and an Aflvenic flux. In conventional flare models, 
such as the CSHKP model, the downward moving newly-closed field lines ultimately join the underlying closed magnetic 
structure. This requires that they slow down, reducing the kinetic energy in the flow. The length of the reconnected flux 
tubes also decreases, and the magnetic energy associated with the current flowing along these field lines also reduces 
as they are dragged into their final location. This implies that both the kinetic energy initially in the outflow, and 
the magnetic energy transported by the outflow decrease. It is assumed here that the deflnition of the energy-release 
region is extended to include this further energy conversion. It follows that the relevant energy emerges as an Alfvenic 
flux, propagating along the newly cl osed fiel d lines. T hese re marks support earlier suggestions that the outflowin g 
energy is Alfvenic (lEmslie &: Sturrock 1982; Haerendell 1994: Fl etcher fc Hudsonll2008l : lKontar. Hannah k Bianll201lD : 
however, there is no detailed model that includes all these effects. A specific model for the Alfvenic energy transport 
away from the energy-release region is presented in the next section. 

4. ALFVENIC ENERGY TRANSPORT 

In this section, a specific (coaxial) model is developed for the Alfvenic transport of energy from the energy-release 
region to the acceleration region. Although highly idealized, this model shows how the large power released can be 
transported over large distances. 

4.1. Setting up the Alfvenic flux 

The magnetic energy inflow into the energy-release region is balanced by an Alfvenic outflow. A complication in 
developing simple models for the inflow and outflow is geometric: the inflow into the energy-release region is modeled 
using cartesian coordinates, and far from the energy-release region the Alfvenic flow is more appropriately modeled 
using cylindrical coordinates. Rather than attempt to formulate a more general geometry that covers both regions, 
cartesian and cylindrical geometries are simply assumed for the inflow and outflow, respectively. 

Consider a cylindrical m odel for a flux tube of radius R with radial coordinate r < R, azimuthal angle </> and axial 
distance z (|Melroselll992| ). The magnetic field is twisted due to the current, and it has components B^,Bz. For 
simplicity the twist is assumed to be small, \B^\ <C \Bz\. The pre-flare current is assumed to be force- free, implying 
Jrt>/ Jz — Brj,/Bz, and the assumption that the twist is small allows one to calculate B^ in terms of Jz and ignore J^. 
It is convenient to write the radial dependence of the current density in the form Jz{r) = </oJi(C): with = r / R a, 
dimensionless radial distance. The function ji(^) describes the current profile in the flux tube prior to the fiare. Prior 
to the flare, there is no electric field, no plasma motion and no energy flux. At the onset of the flare, the inductive 
electric field is set up, and this is identified as a radial field, Er oc g(t). The displacement current, SodEr/dt cx dg{t)/dt, 
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Fig. 2. — A schematic of the coaxial model for Alfvenic energy transport in the upward current region (/q > 0). The dark box on the 
right is the energy-release (ER) region, the central section is the coaxial model, with a current at r < ri , and a current I^Mt = lo ~ hn 
(indicated as a return current) at ri < r < ij; the lightly shaded box on the left is the acceleration region (AR), and on the far left is the 
current, /acc carried by precipitating electrons, and the return current, /ret = lo — lacc from the photosphere. 



drives a radial polarization current, Jr cc dg{t)/dt. This radial current transfers current across field lines, changing the 
current profile as a function of time, concentrating the current density at smaller radial distances, thereby increasing 
i?0 at r < R. 

The radial electric field implies a vortex motion, ~ —E^/B^ oc g{t), where is assumed to be independent 
of time within the cylinder. During the onset of the energy release, increases in magnitude, driven by the force 
— JrBz oc dg{t)/dt. The vortex motion increases the twist, B^jy/Bz- Thus the onset of the energy release launches 
a propagating twist, which may be interpreted as an Alfvenic front that sets up the Alfvenic flux. The front and 
the subsequent energy flux propagate at the MHD speed, vq = va/^^ + v\/cpy^'^. The profile of the front may be 
described by a function g(t — z/vq)^ whose value at z = is df(t)/dt. After the initial transient phase, described by 
dg{t)/dt = d^f{t)/dt^, it is assumed that the energy flux approaches a constant, with g{t) 1. Before discussing the 
Alfvenic flux further, it is appropriate to introduce a specific model for the change in the current profile. 

4.2. Change in axial current profile 

Let the current profile when the Alfvenic flux is established {g{t) ^ 1) be described by j2(C)- As the current profile 
changes at z = 0, it may be described by ~ .9(^)1 + hiOgit)^ with g{t) = prior to the onset of the energy 

release, and with g{t) — 1 after the initial transient phase. At points z > the changing current profile is described 

by 

J. = Mn{m-g{t - z/v^)] +j2{09{t - z/v^)}. (17) 

The axial component of the Poynting vector, which describes the energy flux, is Sz = ErB^/iiQ. The azimuthal 
component of the magnetic field, before and after the Alfvenic fiux is set up, is 



ttR 



t Jo 



(18) 



with Iq — ttR'^Jq. The radial electric field is related to B^ through dEr/dz = dB^/dt, and assuming Er oc g{t — z/vq), 
this gives 

Er=Vo{B^2-B^,)g, (19) 
with g = g{t — z/vq) and -B01.2 = [B^{r)]i^2- The Poynting fiux is then 



Sz 



Mo 



-[5^2 - B^i][gB^2 + (1 - g)B^i]g, 



(20) 



After an initial transient phase, the Alfvenic fiux approaches a constant value vo[Brf,2 — B(j,i]B^2/ fJ-o- 

4.3. Coaxial model 

A simple model for the current profile has ji(^) independent of ^, and j2(C) equal to one constant for ^ < ^1 and 
another constant for ^1 < ^ < 1, as illustrated in Fig. [2] In fact, an arbitrary force- free cylindrical current profile can 
be approximated by a set of concentric rings with the curr ent density a constant within each rin g, with solution for 
B^ involving Bessel functions of the first and second kind ("Me lrose. Nicholls. fc Broderickll 19941 ) . The simple model 
proposed here is a reasonable approximation to such a model for a current profile that is strongly concentrated in an 
inner region, r < Ri — R£,i, provided that the twist is sufficiently small, that is, for B^/B^ <C 1. It is convenient to 
introduce the current, /in, in the inner region, r < and the current, lout, in the outer region, Ri < r < R, with 
I'm + lout = Iq , where Iq is the current in the flux tube prior to the onset of the flare. 

After the onset phase, as g(t) approaches unity, the current density is assumed to approach the constant values a Jq 
for ^ < ^1 and 6J0 for ^1 < ^ < 1. This model involves three parameters, a, &, ^1, only two of which are independent, 
with (a — b)^l = 1 — 6. For g{t) -> 1, this model gives, in terms of (|18l) . 



"'^^^ 1 6e/2 + (l-&)/2^, 



^ <ei, 

< e < 1- 



(21) 
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With bi{() = ^ with (UHl) gives 



Introducing the potential such that = -d^{r)/dr, ([^ with $(0) = implies 



4^ I (i-&)[in(eva + i-r], a<e<i 

The power at r < in the Alfvenic flux is 



2ttr^ dees AO = 

Jo 



Air 

{a-i)e/2{i-b)\ e<a, 
s{o = {, e {i- 2b){e~ei) ng^-et) . (24) 

1-5 2(1-6) ' ^ 

The model involves two independent parameters, b and ^i, with a related to them by (a — b)^f = 1 — 6. There 
is zero Alfvenic flux for 6=1, corresponding to the initial current profile; the power increases with decreasing 6, as 
the current becomes more concentrated near the axis. For 6 = 0, corresponding to all the current in the inner region 
r < ^li?, the power in the Alfvenic flux is of order /o$o with $o = S^/g, which is consistent with the total power 
released in a flare. In this way, the model shows how the released magnetic energy can be transported along field lines 
in a quasi-stationary manner. The parameter 6 can be negative, as discussed in the next section in connection with a 
return current and the number problem. 

The parameter determines the inner radius, ^i-R, and is limited by two requirements. First, the drift speed of the 
electrons, aJo/erie, relative to the ions must not exceed the threshold for the development of an instability that leads 
to the onset of anomalous resistivity. Second, the twist B^p/B^ must not exceed the threshold (of order unity) for the 
current-carrying region to become kink unstable. 

4.4. Inclusion of the displacement current 

The displacement current i s usually ignored in M HD, and it is important to include it in discussing acceleration by 
E\\. This was emphasized bv iSong fc Lvsad (|2006( ). who derived two coupled equations for the t and z derivatives of 
the parallel components of the net current density and the vorticity, O, — curl u. In the cylindrical model assumed 
above, the parallel vorticity is = r~^ d[rurj,{r)\/ dr . 

Neglecting terms that are not important in the present discussion, the equations derived bv iSong fc LvsakI ()2006f ) 
are 

J'.J + ..f, ,25, 

at Ho oz at 

where the net current J' includes the displacement current, and with 

^ v\) dt B, dz- 

On differentiating either (pSj) or I^E^ with respect to t and using the other, one finds that J|' and ilz satisfy the 

equation for an Alfven wave propagating at vq = wa/(1 + ^a/'^^Y^'^ ■ This justifies the assumption made above that 
the temporal variation g(t) at z = propagates away as g{t — z/vq) at z > 0. 

5. ACCELERATION/DISSIPATION REGION 

The acceleration of electrons in a solar flare is an outstanding unsolved problem. In this section, the ingredients 
needed to solve this problem within the framework of the generic model are discussed. The acceleration mechanism is 
assumed to be analogous to the mechanism for auroral electrons, due to -Ey , in a region remote from the energy-release 
region. 

5.1. Requirements on the acceleration region 
The following macroscopic requirements are plausible for the acceleration region. 

1. The Alfvenic flux, incident on top of the region decreases to zero at the bottom of the region. 

2. The potential is across held lines at the top of the region, is zero at the bottom of the region, and is partially 
along field lines within the region, producing . 

3. Electrons are accelerated by to £ = eA$, A$ = $/Af, with M > 1. 
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Fig. 3. — A cartoon of the upward current region of the aurora. T he potential contours (soUd Hnes) indicate a low-altitude and high- 
altitude acceleration region with the auroral cavity in between. [After lErgun et al.l |2004)] 

4. The available power, Iq^q, is transferred to electrons, accelerated at the rate N = MIq/c. 

The most challenging of these constraints is the large value of M, which may be interpreted in terms of the "number 
problem" as the number of times an individual electron would need to be re-accelerated. 

These requirement apply specifically to the upward-current regions above each footpoint, such that the energetic 
electrons propagate downward. As in a substorm, a cceleration also occurs in the nearby downward-current regions 
where electrons are accelerated upward (|Ergun et al.l [2003). 

5.2. Return current and the number problem 

Th number problem (|MacKinnon fc Brownlll989D presents an important difference between flares and substorms. In 
both cases the power is /$ with N = MI /e and e = e^/AI, but with a large multiplicity (M = 10^ with the fiducial 
numbers assumed he re) in flares and M of order unity in substorms. In substorms, inverted- V electrons are accelerated 
in an auroral cavity (jBenson fc Calver3IT9 79'). and as a result of "charge starvation" the flux tube contains essentially 
only the energetic electrons ()Ergun et al.ll200Q.) . In flares, the multiplicity is conventionally interpreted in terms of a 
return current replacing the electrons in the acceleration region M times. An alternative interpretation is that the 
acceleration region and the hard X-ray source are co-located, and the same electrons are accelerated (and lose their 
energy emitting X-rays) a total of M times (Brown et al. 2009). 

If the coaxial model developed above is driven hard enough, it can adjust to maximize the power released. If no 
return current is allowed, 6 > 0, the limiting case is when all the current is confined to the inner region, which becomes 
as small as possible. This corresponds to Jin = Iq, lout — 0, a^^ = 1, & = 0. When a return current is allowed, b < 0, 
an even larger power is possible. For 6 < 0, /out has the opposite sign to Im, and corresponds to a return-current for 
^ > ^1, as illustrated in Fig. [2] It follows from (|24|) that the power increases with increasing |&| for 6 < 0. The model 
becomes untenable for b large and negative. Nevertheless, this case gives an indication of how the number problem 
can be resolved. To show this, let us assume b ~ —M, and explore the implications. 

According to (P^ . the power in the Alfvenic flux increases oc (l — b)^ « AP ioi —b ~ M ^ 1. This can be attributed 
to the combination of two effects. First, the total current flowing in the region r < ^iR is much greater than Iq, with 
(all but Jo of) this direct current being provided by recycling the return current, which closes across field lines in the 
photosphere and in the energy-release region. Second, the large current /in ~ MIq implies a correspondingly large 
cross-field $ « Af$o, which is implicit in the assumption of Alfvenic fiux. With M w 10^, the extremely large values 
of / and 3" are not tenable. In the coaxial model the direct current flows in a single current channel, at r < ^iR, and 
relaxing this assumption opens the possibility of formulating a model that can resolve the number problem. 

Suppose that the direct current is of order MIq flowing in M parallel channels, each carrying a current Iq. The 
return current flows between these channels, and is of magnitude (M — l)/o- If the potential associated with each 
channel is $0/-^^, with $0 = ^aIq, the Alfvenic flux in each channel transports a power Iq^q/AI, giving the total 
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power /o$o after summing over the M channels. In the acceleration region, this power is transferred to electrons with 
energy e = e^o/M at a rate Iq/c in each channel, resulting in MIo/e electrons per second with this energy from the 
M channels. The implied cross-field potential in any given channel is ^q/M w 10^ V, consistent with the observed 
energy, e, of order lOkeV. 

Although these arguments are heuristic, they give a clear indication of how the number problem can be resolved. A 
detailed model incorporating these ideas is needed to make further progress. 

5.3. Development of 

An understanding of acceleration by involves two complementary aspects: the macroscopic electrodynamics 
associated with the development of E'n , and the microphysics associated with the specific (local and transient) structures 
with fi'ii ^ 0. Acceleration by is widely accepted for auroral electrons in a substorm (i Block fc FalthammaiifTOQOl) . 
where E\\ can be observed in situ (jMozer fc Kle tzing 1998), and also for lo- related electrons in Jupiter's magnetosphere 
(|Su et al.l 120031 ). However, the microphysics involved in setting up and maintaini ng E\\ is inadequately understood, 
with none of the many suggested mechanisms being entirely satisfactory (|Borovskvlll993D . As remarked in paper 1, in 
auroral acceleratio n E\\ is highly structured on the microscale, with individual structures being modeled as, for example, 
phase-space holes (|Ergun et al.lll998bt ) and kinetic Alfven waves (jChaston et al.|[2003D . Discussion of acceleration by E'n 
in the solar corona also emphasizes microphysics (|Bian. Kontar. fc Brown|[2010[ ). The viewpoint adopted here is that 
important aspects of acceleration by E^^ can be understood in terms of the macroscopic electrodynamics, independent 
of the microphysics that sustains This view was expressed bv iHaerendell (|1994f) . who argued that the appeal to 
microphysics (double layers and anomalous resistivity) is required only to support other arguments that require the 
appearance of . 

One argument that an E\\ must appear is that the downward Alfvenic fiux, transporting the released magnetic energy, 
has a large Ex and associated vortex motion, u, and these cannot be sustained through the denser chromospheric 
and photospheric plasma. As illustrated in Fig. [3] for the auroral application, the surfaces of constant potential, 
corresponding to Ej_, must close above the conducting ionosphere or chromosphere, implying that must develop. 
Another argument is that line-tying, associated with the inertia of the denser regions of the ionosphere or chromosphere, 
precludes the twisting motion driven from above being maintained. The energy in the Alfvenic flux must be reflected, 
absorbed or converted into another form somewhere above the boundary where line-tying becomes effective. This 
implies that along a given held line, ^ must reduce to zero, and this violates the frozen-in condition. The 
development of E\\ is ess ential to allow a decoupling between two regions in relative motion perpendicular to B, called 
fracturing bv IHaerendell (|1994 [2007[ l. These arguments imply that there must be a region with E\\ =/= irrespective of 
the speciflc microphysics that sustains . In the model suggested above in connection with the return current and the 
number problem, the cross-fleld potential is ^q/M, and the implication is that this leads to a parallel potential ^o/M 
along held lines in the acceleration region. The mocrophysics involved in how E^^ develops is not directly relevant: the 
system must find some form of microphysics to satisfy the macroscopic requirements. 

6. DISCUSSION AND CONCLUSIONS 

In this paper, and in paper 1, a generic model for magnetic explosions is formulated, based on magnetospheric 
substorms, and applied to solar fiares. Three ideas are central to the generic model: the separation into different scales 
(global, macro, micro) and into different (energy-release, acceleration, current-closure) regions; the insistence that the 
intrinsic time dependence be included explicitly (through Ejnd, Ujnd, Jdisp, Jpoi); and, that the electron acceleration 
is due to a parallel electric field, E\\ . 

What triggers a flare is an unsolved problem. The recognition that the inflow into the energy-release region is driven 
by the changing magnetic field offers a new way of thinking about the trigger. Prior to a flare the magnetic structure 
is in a steady state, with any pressure gradients balanced by the J x B force, with a very small angle between J and 
B. Once a flare is triggered, the unbalanced Maxwell stress, mediated by Ejnd and Jdisp, drives the energy release, 
without requiring a mechanical force. The trigger can be a perturbing pressure gradient that initiates inflow into 
the energy-release region, causing B to start to change with time. The changing B can then takes over, driving the 
magnetic energy into the energy-release region and powering the explosion. 

It is assumed that the energy release occu rs in two stages, involving two widely separated regions with an Alfvenic 
flux transporting the energy between them (|Fletcher fc Hudsonll2008f ). The energy transport is modeled in terms of 
an Alfvenic energy flux associated in a current-carrying flux tube. A coaxial model is developed, in which the direct 
current becomes strongly concentrated near the axis of the cylinder. A limiting case of this model is when all the 
current is confined to narrow channel, and in this limit the total power released can be transported Alfvenically from 
the energy-release region to the acceleration region. 

The "number problem" has presented a formidable challenge to all theories of electron acceleration in flares: the 
hard X-ray date imply either that the number of electrons accelerated is M times the number available, or that the 
available electrons are re-accelerated M times, with M of order 10^. The observations required electrons with an 
energy e^o/M accelerated at a rate N = MIo/e, such that the total power, Iq^o, is that released in the flare. Based 
on the coaxial model developed here, it is suggested that the number problem may be resolved by incorporating three 
changes. First, the single direct current channel is replacement by M current channels, each carrying a current Iq. 
Second, there is a return current of magnitude (M — l)/o that flows between these channels, with the current lines 



Inductive electric field 



9 



closing in the photosphere and the energy-release region. Third, the cross-field potential associated with each channel 
is ^q/M^ and this potential develops along field lines in the acceleration region. 

Another feature of the model is that there are both upward and downward current regions at each footpoint of the 
flaring flux tubes. The energetic electrons that produce hard X-rays as they precipitate are attributed to acceleration in 
the upward current region, and the electrons that escape to produce type III radio bursts are attributed to acceleration 
in the downward current region. These electrons do not come from a single distribution and need not have the same 
properties. Energetic ions, accelerated downward in the downward current region, arc necessarily slightly displaced 
from the electrons accelerated downward in the upward current region. Such effects need to be treated quantitatively 
in a detailed model. 
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